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A Radiation  Energy-Angle  Algorithm 
for  Use  in  Personnel  Dosimetry 

Margarete  Ehrlich 

National  Institute  of  Standards  and  Technology 
GaUhersbiirg,  MD  20899 


In  an  earlier  publication,  the  author  described  an  algorithm  for  arriving  at  the  value  of 
the  individual  dose  equivalent  from  the  readings  of  dosimeters  having  two  areas  of  different 
energy-  and  angle-response  functions.  The  algorithm  is  briefly  reviewed  here,  and  suggestions 
are  made  for  obtaining,  by  numerical  means,  a plausible  average  of  all  solutions  possible  for  a 
given  type  of  dosimeter  and  number  of  calibration  data. 

Key  words:  algorithm;  angular  dependence;  energy  dependence;  individual  dose  equivalent; 
personnel  dosimetry;  plausible  average  of  solutions;  response  functions;  response  ratios;  sur- 
faces in  three  dimensions;  three-dimensional  spline  fits. 


1.  Introduction 

In  the  days  prior  to  the  existence  of  the  “effective  dose  equivalent”  and  the  “individual 
dose  equivalent,”  personnel  dosimeters  were  constructed  with  only  little  attention  to  the 
angular  dependence  of  their  response.  Dosimeter  calibration  was  performed  as  if  the 
dosimeters  were  to  be  used  solely  under  perpendicular  incidence  of  both  the  primary  and 
scattered  radiation,  and  dosimeter  evaluation  was  in  terms  of  quantities  that  ignored  the 
asymmetries  in  the  absorption  of  radiation  in  the  human  body.  This  situation  changed  with 
the  advent  of  the  effective  dose  equivalent  and  of  some  of  the  new  operational  dose-equiv- 
alent quantities  [1],  which  take  these  differences  into  account.  As  a consequence,  more 
studies  are  being  performed  on  angular  dependence  of  the  response  of  personnel  dosimeters 
in  current  use.  (See,  e.g.,  Marshall  et  al.  [2],  King  [3]  and  the  various  papers  on  this  topic 
presented  at  two  seminars  on  radiation  protection  quantities  for  external  exposure  held  in 
1985  [4]  and  1988  [5].)  Also,  the  revised  personnel-dosimetry  performance  test  document, 
ANSI  N-13.11  [6],  will  go  far  beyond  the  old  edition  [7],  which  recommended  only  that 
results  of  such  studies  be  available. 

While  new  dosimeter  construction  probably  will  take  angular-dependence  require- 
ments into  account,  in  many  instances  the  use  of  older  types  of  dosimeters  probably  will 
continue  for  some  time.  This  report  deals  with  the  evaluation  of  readings  on  some  of  these 
older  types  of  dosimeters  in  terms  of  either  air  kerma  or  dose  equivalent  by  a procedure  that 
may  make  the  replacement  of  these  dosimeters  altogether  unnecessary.  A brief  review  is 
given  of  this  evaluation  procedure  (“evaluation  algorithm”)  which  was  first  presented  in 
1988  at  the  Seminar  on  Radiation  Protection  Quantities  mentioned  above  [8],  and  steps  are 
suggested  for  improving  it. 
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2.  History 


In  the  early  days,  photographic  emulsions  were  essentially  the  only  radiation-sensitive 
elements  employed  for  personnel  dosimetry  — with  personnel  dosimeters,  in  general,  incor- 
porating both  nuclear-track  emulsions  for  neutron  dosimetry  and  conventional  emulsions 
for  photon  and  beta-particle  dosimetry.  The  response  functions  of  photographic  emulsions 
are  known  to  be  strongly  dependent  on  radiation  energy  and  angle  of  radiation  incidence.^ 
Nevertheless,  the  customary  way  to  interpret  response  to  photons  or  beta  particles  in  terms 
of  a quantity  related  to  the  radiation  exposure  on  an  individual  was  to  assume  perpendicular 
radiation  incidence  upon  the  front  face  of  the  dosimeter.  After  isolating  the  effect  of  neu- 
trons, discrimination  between  the  effects  of  beta  particles  and  photons  of  different  energies 
was  accomplished  by  incorporating  into  each  personnel  dosimeter  a series  of  different 
filters  or  filter  combinations,  covering  most  of  the  surface  area  of  the  conventional  emul- 
sions, and  providing  one  or  two  areas  of  different  response  to  beta  particles  and  photons, 
and  at  least  two  further  areas  of  different  response  to  photons  of  different  energies.  (See, 
e.g.,  the  review  given  in  IAEA  Safety  Series  No.  8.  [9])  The  resulting  evaluation  procedure 
used  for  photon  or  beta-particle  dosimetry  is  described  below,  with  contemporary  nomen- 
clattire. 

2.1  Dosimeter  Calibration  for  Photon  Dosimetry  (Historical) 

The  dosimeters  were  calibrated  solely  with  radiation  of  perpendictilar  incidence.  Cali- 
bration involved  the  following  steps: 

(a)  irradiating  the  dosimeters  in  known  radiation  fields  covering  the  entire  range  of 
interest  of  either  monoenergetic  photons  of  energy,  E,  or  of  narrow  spectral  bands  of 
photons  of  an  average  energy,  E,  at  known  levels  of  the  dosimetric  quantity,  in  terms 
of  which  dosimeter  readings  were  to  be  interpreted; 

(b)  measuring  optical  density,  M,  in  all  emulsion  areas  of  all  dosimeters,  and  plotting 
or  tabulating  the  values  of  M as  a function  of  H,  with  E as  the  parameter. 

(c)  determining  the  response  functions,  M,(E)/if(E),  under  the  ith  filter,  similarly 
Mj(E)IH{E)  under  the ;th  filter,  and  so  on  for  all  filter  combinations  and  all  energies,  E; 
then  taking  the  ratios,  M,(E)/A//(E)  etc.,  of  these  response  functions,  and  plotting  or  tabu- 
lating them  as  a function  of  E for  ail  calibration  energies. 

2.2  Algorithm  for  Evaluating  H from  Measurements  on  the  Field  Dosimeters  Irradiated  with 
Protons  (Historical) 

Evaluation  of  H on  the  field  dosimeters  involved: 

(a)  measuring  the  values  of  Mi,  Mj  etc.,  under  each  filter  of  each  field  dosimeter,  and 
forming  the  response  ratios,  M,  /My  etc.; 

(b)  with  the  aid  of  the  calibration  relationships  obtained  for  the  response  ratios  in  step 

(c)  of  section  2.1,  finding  the  photon  energy  (or,  in  general,  energies),  E,  corresponding  to 
each  of  these  ratios; 


‘ “Response  function”  is  the  customary  term  for  “response  as  a function  of  radiation  energy,”  where  “response” 
is  defined  as  dosimeter  reading  per  unit  of  the  chosen  dosimetric  quantity. 

* In  this  report,  the  dosimetric  quantity  is  represented  by  the  symbol,  H,  regardless  of  whether  it  is  the  old-time 
“ejposure,”  air  kerma,  or  one  of  the  operational  dose-equivalent  quantities. 
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(c)  with  the  aid  of  the  calibration  relationships  between  M and  H obtained  in  step  (b) 
of  section  2.1  with  the  parameter  finding  the  values  of  H corresponding  to  each  mea- 
sured value  of  Mi,  Mj  etc.,  for  the  energies  determined  in  the  preceding  step. 

This  algorithm  seems  more  straightforward  than  it  really  is,  because,  in  general,  both 
step  (b)  and  the  fact  that  the  energy  spectrum  encountered  in  the  field  does  not  match  the 
one  used  for  the  calibration,  may  result  in  more  than  one  energy  value.  The  final  selection 
of  a value  for  H therefore  necessitates  a judicious  choice. 

In  the  following,  it  will  be  shown  how  the  algorithm  described  above  can  be  general- 
ized to  apply  to  the  current  problem  of  determining  the  value  of  the  chosen  quantity,  say 
the  individual  H,  with  the  angle  of  radiation  incidence,  0,  being  assumed  to  vary  between 
at  least  0”  and  90°. 


3.  Generalization  of  the  Historical  Method 

This  generalization  should  be  applicable  to  all  types  of  radiation  for  which  one  uses 
personnel  dosimeters  that  (a)  have  response  functions  depending  strongly  on  radiation  en- 
ergy and  angle  of  radiation  incidence,  and  (b)  incorporate  at  least  two  areas  of  different 
energy-  and  angle-response  functions.  It  should  be  applicable  also  to  instances  in  which  two 
types  of  detectors  with  different  response  functions  are  used  either  in  the  same  dosimeter  or 
in  two  dosimeters  worn  side  by  side.  Although  not  explicitly  treated,  this  case  is  implied  in 
all  the  following  procedural  steps,  in  which  reference  is  made  to  “detectors  / and The 
method,  which  does  not  require  separate  knowledge  of  the  energy  spectrum  and  the  angu- 
lar distribution  of  the  radiation  encountered  in  the  field,  was  described  by  the  author  in  a 
presentation  made  at  the  1988  Radiation  Protection  Dosimetry  Seminar  [8].  In  this  presenta- 
tion, relatively  encouraging  results  were  shown  of  the  author’s  attempt  to  apply  the  al- 
gorithm in  a rough  and  unsophisticated  manner  to  one  particular  dosimeter  type  in  current 
use.  At  the  same  Seminar,  a similar  idea  for  an  algorithm  was  touched  upon  also  in  a 
presentation  from  the  Physikalisch-Technische  Bundesanstalt  [10]. 

In  the  following,  a brief  review  is  given  of  the  quantities  employed,  and  of  the  calibra- 
tion method  and  the  evaluation  algorithm.  This  is  followed  by  suggestions  of  ways  to 
reduce  the  imcertainties  arising  from  the  multiplicity  of  the  results  obtained. 

3.1  Quantities 

Here,  the  dosimeter  indication  M applies  to  any  type  of  radiation  effect  measured  for 
any  type  of  dosimeter.  The  quantity  H is  taken  to  be  the  individual  dose  equivalent  (shallow 
or  penetrating)  [1],  which  is  still  awaiting  its  final  definition  [11,12].  Both  M and  H are  taken 
to  be  functions  of  radiation  energy,  £,  and  of  angle  of  radiation  incidence,  0,  with  H{E,%) 
being  known  at  present  at  least  approximately  from  computations  of  the  directional  dose 
equivalent.  (The  quantity  H could  have  been  taken  also  to  represent  air  kerma,  and  then 
could  have  been  designated  by  H{E).) 

3.2  Dosimeter  Calibration  (Generalized) 

For  dosimeters  incorporating  two  detectors,  i and  j,  with  different  response  functions, 
dosimeter  calibration  ent^s: 

(a)  irradiating  the  dosimeters  in  known  radiation  fields  of  the  type  of  radiation  under 
consideration,  covering  the  entire  energy  and  angle  range  of  interest  and  known  levels  of 
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the  quantity,  H; 

(b)  measuring  dosimeter  indication,  M,  as  a function  of  H,  for  detectors  i and  j,  over 
the  range  of  energies  of  interest  and  over  the  range  of  angles  of  radiation  incidence  between 
at  least  0°  and  90®,  and  tabulating  or  plotting  the  values  of  M as  a function  of  H for  all 
values  of  E and  0; 

(c)  forming  the  response  functions,  Mi{E,Q)IH(E,Q)  and  Mj{E,Q)/H{E,Q),  and  tabu- 
lating or  plotting  them  against  both  E and  0 for  all  values  oli,  j,  E and  0;  and 

(d)  forming  the  response  ratios, 

{Mi{E,Q)IH{E,Q)}l{Mj{E,Q)IH{E,^)}, 

and  tabulating  or  plotting  them  against  both  E and  0 for  all  values  oii,  j,  E and  0. 

The  data  may  be  plotted  either  in  two  dimensions,  resulting  in  families  of  curves  with 
either  £ or  0 as  the  parameter,  or  m three  dimensions,  resulting  in  a family  of  surfaces. 
Several  types  of  numerical  fitting  procedures,  e.g.,  polynomial  least-squares  fits  or  various 
types  of  spline  fits,  are  available  for  use  in  either  case. 

Examples  of  families  of  curves  in  two  dimensions,  for  only  four  energies  and  three 
angles,  drawn  simply  by  eye,  were  included  in  the  earlier  presentation  [8].  For  some  of  the 
same  data,  examples  of  surfaces  in  three  dimensions  for  M,  IH  versus  £ and  0 and  for  Mi  IMj 
versus  E and  0 were  produced  by  Boisvert  [13]  with  a computer  program  on  the  NIST 
mainframe  computer.  These  surfaces  are  shown  in  figures  1,  2 and  3. 
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Figure  1.  Example  of  three-dimensional  spline  of  a response  function.  Mi  IH  [13]. 
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Figure  2.  Example  of  three-dimensional  spline  of  a response  function,  Mj/H  [13]. 
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Figure  3.  Example  of  three-dimensional  spline  of  a response-function  ratio,  MilMj  [13]. 
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33  Algorithm  for  Evaluating  H from  Measurements  on  the  Field  Dosimeters  (Generalized) 

Evaluation  of  H on  the  field  dosimeters  entails: 

(a)  measuring  values  of  M for  the  detectors  i and  j of  each  field  dosimeter,  and 
forming  the  response  ratios,  Mi/Mj; 

(b)  finding  the  (E  ,0)  pair  (or,  in  general,  pairs)  that  correspond  to  the  values  of  Mi  IMj 
obtained  in  step  (d)  of  section  3.2,  on  the  calibration  surfaces  or  curves  representing 
Mi{E ,Q)IMj{E ,Q)  as  a function  of  E and  0; 

(c)  obtaining  the  values  (or  sets  of  values)  of  MilH  and  MjlH  for  each  set  of  (E,0) 
pairs  on  all  field  dosimeters,  with  the  aid  of  the  calibration  data  of  M,(E,0)/H  and  M;(E,0)/ 
H,  respectively  (step  (c),  sec.  3.2);  and 

(d)  determining  the  value  of  H corresponding  to  each  measured  value  of  Mi  IMj  from 
the  calibration  relationships  between  Mi,  Mj  and  H (step  (b)  of  sec.  3.2). 

3.4  Multiplicity  of  Results 

Step  (b)  of  the  algorithm  outlined  in  section  3.3  generally  results  in  a set  of  different 
(£,0)  pairs,  all  corresponding  to  the  same  measured M, /My  ratio.  Their  number  is  likely  to 
increase  with  the  number  of  calibration  energies  and  angles;  but  so  is  the  Ukelihood  that 
they,  as  a set,  reflect  more  accurately  the  actual  field-irradiation  conditions— for  which 
energy  spectra  and  angular  distributions  may  be  very  different  from  those  used  during 
calibration.  Also,  it  is  to  be  expected  that  the  set  of  the  more  reliable  results  will  stem  from 
values  of  M measured  on  the  detector  i or  j for  which  the  response  functions  vary  more 
strongly  with  E and  0. 

3.5  Assigning  Weights  to  the  Resulting  Individual  H Values 

One  then  needs  to  decide  on  a method  for  arriving  at  a unique  value  of  H representing 
the  set  of  H values  obtained  in  step  (d).  In  first  approximation,  one  may  consider  taking  the 
algebraic  (unit-weight)  mean  of  the  H values  in  each  set.  This  is  the  simple-minded  ap- 
proach for  which  results  for  the  more  E - and  0-dependent  detector  were  shown  in  the 
original  presentation.  [8]  Following  is  a detailed  discussion  of  the  mechanics  involved  in 
some  of  the  individual  steps  of  the  proposed  algorithm,  suggesting  a more  plausible  averag- 
ing method.  The  basic  idea  for  this  method,  which  as  yet  has  not  been  tried,  originated  with 
M.  Danos  [14],  and  was  later  elaborated  on  in  discussions  with  S.  M.  Seltzer  [15]. 

3.5.1  Mechanics  Involved  in  Arriving  at  (E,0)  Values 

Let  us  assume  that  the  calibration  surfaces  representing  Mi{E,Q)IH{E,Q),  My(E,0)/ 
H(£,0),  and  M/(E,0)/My(E,0)  as  functions  of  E and  0 (sec.  3.2)  have  been  fitted  numeri- 
cally, and  that  coordinates  of  surface  points  for  equidistant  values  along  the  E and  0 axes 
(“grid  points”)  have  been  obtained.  Let  us  further  assume  that  the  response  ratios.  Mi  IMj, 
have  been  determined  on  the  field  dosimeters  by  step  (a)  of  the  algorithm  (sec.  3.3).  One 
then  is  ready  for  step  (b),  which  entails  cutting  the  calibration  surfaces  representing 
Mi{E,Q)IMj{E,Q)  as  a function  of  £ and  0 by  planes  given  by 


Mi  IMj  = constant. 
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where  M,-  IMj  is  one  of  the  ratios  derived  from  the  measurements.  This  step  results  in  a 
curve  of  intersection,  with  points  representing  the  set  of  (£,0)  pairs  for  the  particular 
Mi  IMj  ratio.  The  subsequent  steps  (c)  and  (d)  then  may  be  carried  out  for  all  the  points 
along  the  length,  /,  of  the  curve  of  intersection,  resulting  in  a set  of  H values  corresponding 
to  these  points. 

3.5.2  Averaging  the  H Values  Along  the  Curve  of  Intersection 

Let  us  pick  a set  of  (£,0)  pairs  along  the  curve  of  intersection  of  the  plane, 
Mj /My  — constant  with  the  surface  representing  Mi(E,0)/My(E,0)  as  a function  of  .E  and  0, 
that  has  equidistant  spacing  of  either  the  E or  the  0 values,  i.e.,  that  is  either  of  the  form 


or  of  the  form 


(E,0)4=(E,0i), 


where  it  is  a running  index.  Each  of  these  points  can  be  associated  with  a length.  A/*,  along 
the  curve,  where  A/^  may  be  approximated  as 


Ai*  ~ ^ {Ek+i—  EkY + {Qk+i~  • 

Let  us  then  designate  by  H(Ek,Qk)  the  value  of  H obtained  in  step  (d)  for  the  length 
increment.  A/*,  of  the  curve  of  intersection.  A plausible  value,  representing  the  set  of 
H{Ek,Qk)  values  along  the  line  /,  would  be  the  weighted  average  given  by 

Hp- XH(E„0OA/*/iM*. 


Alternatively,  one  might  try  to  obtain  similarly  weighted  averages  of  the  values  of  E*  and 
0jt  along  the  curve  of  intersection  obtained  in  step  (b),  and  then  proceed  with  steps  (c)  and 
(d)  only  for  the  one  (Ep,0p)  pair. 


4.  Accuracy  of  the  Method 

At  this  stage,  it  is  difficult  to  predict  how  accurate  the  results  will  be  that  one  can 
obtain  with  the  described  algorithm.  Among  other  things,  accuracy  will  depend  on  such 
factors  as  the  shape  of  the  response  functions  and  response-function  ratios  of  a particular 
type  of  dosimeter,  and  how  well  these  functions  are  known— i.e.,  on  how  many  calibration 
points  are  available.  Once  the  numerical  procedures  are  worked  out,  the  algorithm  should 
be  tested  on  a type  of  dosimeter  that  fulfflls  the  criteria  given  in  section  3,  for  which  data 
are  available  for  a sufficient  number  of  radiation  energies  and  angles  of  radiation  incidence. 

The  test  could  be  carried  out  with  photons,  neutrons  or  beta  particles,  or  a mixture  of 
them.  In  the  latter  case,  the  assumption  could  be  made  that  the  discrimination  between  the 
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types  of  radiation  has  been  carried  out  successfully  beforehand.  Also,  initially,  the  test 
could  be  run  as  if  the  calibration  dosimeters  were  identical  with  the  field  dosimeters  for 
which  irradiation  levels,  energies  and  angles  are  unknown.  But  this  leads  to  an  inherent 
advantage  regarding  the  compatibility  of  calibration  and  field  conditions.  Therefore,  the 
test  must  be  extended  also  to  “field  dosimeters”  irradiated  at  energies  and  angles  different 
from  those  used  during  calibration. 
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